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Abstract 

Given a finite set of arbitrarily distributed points in affine space with arbitrary multiphcity structures, 
we present an algorithm to compute the reduced Grobner basis of the vanishing ideal under the lexico- 
graphic ordering. Our method discloses the essential geometric connection between the relative position of 
the points with multiplicity structures and the quotient basis of the vanishing ideal, so we will explicitly 
know the set of leading terms of elements of /. We split the problem into several smaller ones which can 
be solved by induction over variables and then use our new algorithm for intersection of ideals to compute 
the result of the original problem. The new algorithm for intersection of ideals is mainly based on the 
Extended Euclidean Algorithm. 

Keywords: vanishing ideal, points with multiplicity structures, reduced Grobner basis, intersection of 
ideals. 



1. Introduction 

To describe the problem, first we give the definitions below. 

Definition 1: D C Nq is called a lower set as long as Vd G D if 7^ 0, d — lies in D where 
= (0, . . . , 0, 1, 0, . . . , 0) with the 1 situated at the i-th position (1 < i < n). For a lower set D, we define 
its limiting set E{D) to be the set of all € Ng — D such that whenever /3i 7^ 0, then (3 — Ci G D. 

As showed in Fig.l below, there are three lower sets and their limiting sets. The elements of the lower 
sets are marked by solid circles and the elements of the limiting sets are marked by blank circles. 
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Fig.l: Illustration of three lower sets and their limiting sets. 

Let k he a field and p be a point in the affine space fc", i.e. p — {pi, . . . ,pn) G fc". Let k[X] be the 
polynomial ring over k, where we write X — (Xi, X2, . . . , Xn) for brevity's sake. 

Definition 2: {p, D) represents a point p with multiplicity structure D, where p is a point in affine 
space fc" and Z? is a lower set. (jZ? is called the multiplicity of point p (here we use the definition in [3]). 
For each d — (di, . . . ,dn) G D, we define a corresponding functional 

adi + ...+d„ 

Hence for any given finite set of points with multiplicity structures H = {{pi, Di), . . . , {pt,Dt)}, we 
can define m functionals where m = (jZ3i + . . . + tt_Df. Our aim is to find the reduced Grobner basis of 
the vanishing ideal 1(H) = {/ e k[X];Li{f) — Q,i — l,...,m} under the lexicographic ordering with 
Xi ^ X2 ^ . . . ^ X„. 

There exists an algorithm that provides a complete solution to this problem in [4]. However, our answer 
for the special case of lexicographical ordering will be in a way more transparent than the one above. The 
ideas are summed-up as follows: 
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• Construct the reduced Grobner basis of I{H) and get the quotient basis D{H) by induction over 
variables. 

• Get the quotient basis D{H) purely according to the geometric distribution of the points with 
multiplicity structures. 

• Split the original problem into smaller ones which can be converted into 1 dimension lower problems 
and hence can be solved by induction over variables. 

• Compute the intersection of the ideals of the smaller problems by using Extended Euclidean Algo- 
rithm. 

There are several publications which have a strong connection to the work presented here. Paper [5] 
give a computationally efficient algorithm to get the quotient basis of the vanishing ideal over a set of 
points with no multiplicity structures and the authors introduce the lex game to describe the problem. 
Paper [6] offers a purely combinatorial algorithm to obtain the linear basis of the quotient algebra which 
can handle the set of points with multiplicity structures but it does not give the Grobner basis. For a finite 
set of points with multiplicity structures, our algorithm obtains a lower set by induction over variables and 
constructs the reduced Grobner bases at the same time. It is only by constructing Grobner basis we can 
prove that the lower set is the quotient basis. 

One important feature of our method is the clear geometric interpretation, so in Section 2 an example 
together with some auxiliary pictures will be given in the first place to demonstrate this kind of feature 
which can make the algorithms and conclusions in this paper easier understood for us. In Section 3 and 4, 
some definitions and notions are given. Section 5 and 6 are devoted to our main algorithms of computing 
the reduced Grobner basis and the quotient basis together with the proofs. In Section 7 we demonstrate 
the algorithm to compute the intersection of two ideals and some applications. 



2. Example 



First we give two different forms to represent the set of points H with multiplicity structures. 

For easier description, we introduce the matrix form which consists of two matrices {V = {Pi,j)mxn,1^ = 
(rfij)mxn) with VijVi denoting the i-th row vectors of V and V respectively. Each pair {Vi,'Di} (1 < i < 
m) defines a functional in the following way. 



Li{f) 



,l+...+di. 



di, 



Xl=Pi,l,...,X„=Pi 



OX^ ' ... OXn 

And the functional set defined above is the same with that defined by H in Section 1. 

For example, given a set of three points with their multiplicity structures {{p-i,Di), {p2, i?2), (ps, £^3)}, 
where p, = (1, l),p2 = (2, 1),P3 = (0, 2), D, = {(0, 0), (0, 1), (1, 0)}, D2 = {(0, 0), (0, 1), (1, 0), (1, 1)}, D3 = 
{(0, 0), (1, 0)}, the matrix form is like the follows. 



/ 1 1 \ 

1 1 

1 1 

2 1 



r = 



\0 2J 



/ \ 

1 

1 



1 

1 

1 1 



V 1 y 



For intuition's sake, we also represent the points with multiplicity structures in a more intuitive way 
as showed in the left picture of Fig. 2 where each lower set which represents the multiplicity structure of 
the corresponding point p is also put in the afHne space with the zero element (0,0) situated at p. This 
intuitive representing form is the basis of the geometric interpretation of our algorithm. 

We take the example above to show how our method works and what the geometric interpretation of 
our algorithm is like: 

Step 1: Define mapping tt : H 1-^ k such that (p — (pi, . . . ,pn),D) ^ H is mapped to p„ € k. So 
H = {(pi, Di), (p2, -D2), (P3, D^)} consists of two 7r-fibres: = {(pi, Di), (p2, -D2)} and H2 = {(ps, D^)} 
as showed in the middle and the right pictures in Fig.2. Each fibre defines a new problem, so we split the 
original problem defined by H into two small ones defined by Hi and H2 respectively. 
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Fig. 2: The left picture represents H. The middle one is for Hi and the right one for H2- 

Step 2: Solve the small problems. Take the problem defined by Hi for example. 
First, it's easy to write down one element of I{Hi): 

h = {X2 - 1){X2 - 1) = {X2 - If e I{Hi). 

The geometry interpretation is: we draw two lines sharing the same equation of X2 — 1 = to cover 
all the points as illustrated in the left picture in Fig. 3 and the corresponding polynomial is /i. 



Fig. 3: Three ways to draw lines to cover the points. 

According to the middle and the right pictures in Fig. 2, we can write down another two polynomials 
in I{Hi): 

/2 = {X2 - l){Xi - l){Xi - if and /3 = {Xi - \f{Xi - 2f . 

It can be checked that Gi = {fi, /2, /s} is the reduced Grobner basis of I{Hi), and the quotient basis is 
{1, Xi,X2, XiX2,Xi, X2X1, Xf}. In the following, we don't distinguish explicitly an n-variable monomial 
j^di^da _ ^ jj-d„ ^-^-^Yi the element (di, ^2, ■ • • , dn) in Nq , and we denote the quotient basis of I{H) by D{H). 
Hence D{Hi) can be written as a subset of NJf : {(0, 0), (1, 0), (0, 1), (1, 1), (2, 0), (2, 1), (3, 0)}, i.e. a lower 
set, denoted by D . 

In fact we can get the lower set in a more direct way by pushing the points with multiplicity structures 
leftward which is illustrated in the picture below (lower set D is positioned in the right part of the picture 
with the (0,0) element situated at point (0,1)). The elements of the lower set D in the right picture in 
Fig. 4) are marked by solid circles. The blank circles constitute the limiting set E{D ) and they are the 
leading terms of the reduced Grobner basis {/i, /2, /a}- 



1 -i> O 00 O 



1 2 Xi 

Fig. 4: Push the points leftward to get a lower set. 



In the same way, we can get the Grobner basis G2 = {/ii,^2} and a lower set D for the problem 
defined by H2, where = (X2 - 2),/i2 = Xf,D" = {(0, 0), (1, 0)}. 

Step 3: Compute the intersection of the ideals I{Hi) and I(H2) to get the result for the problem 
defined by H. 

First, we construct a new lower set D based on D , D in an intuitive way: let the solid circles fall 
down and the elements of D rest on the elements of D to form a new lower set D which is showed in the 
right part of Fig. 5 and the blank circles represent the elements of the limiting set E{D). 
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Fig. 5: Get the lower set D based on D and D 



Then we need to find '^E[D) polynomials vanishing on H with leading terms being the elements of 
E{D). Take XIX2 G E{D) for example to show the general way we do it. 

We need two polynomials which vanish on Hi and H2 respectively, and their leading terms both have 
the same degrees of Xi with that of the desired monomial X1X2 and both have the minimal degrees of 
X2. It's easy to notice that /2 and Xi ■ /12 satisfy the requirement and then we multiply /2 and Xi ■ /12 
with hi , /i respectively which are all univariate polynomials of X2 to get two polynomials qi , q2 which 
both vanish on H. 

91 - /2 • /ii - {X2 - l){Xi - l){Xi - 2f{X2 - 2), 



92 



Xi-h2- h^Xl{X2-l) 



Next try to find two univariate polynomials of X2'. ri , r2 such that gi • + (72 • r2 vanishes on H (which 
is apparently true already) and has the desired leading term XiX2- 

To settle the leading term issue, write (71,92 as univariate polynomials of Xi. qi = {X2 — 2)(X2 — 
l)Xl - (5X| - 15X2 + 10)Xf + (8X| - 24X2 + 16)Xi - 4X| + 12X2 - 8,(72 = {X2 - ifXl. Because 
X2 -< Xi and the highest degrees of Xi of the leading terms of 91, (72 are both 3, we know that as long as 



the leading term of {X2 
is also XiX2. 



2)(X2-l)X3.ri + (X2 



1)2X3 . 



r2 is XfX2, the leading term of qi ■ ri + 92 • ^2 



(X2 - 2){X2 - i)xf ■ n + {X2 - ifxf ■ ,'2 

= XUX2 - 1) ((X2 - 2) • n + {X2 - 1) • r2) 

Obviously if and only if (X2 — 2) • ri + {X2 — 1) • ?'2 = 1 we can keep the leading term of 91 • ri + (72 • r2 to 
be XiX2- In this case ri = —1 and ^2 = 1 will be just perfect. In our algorithm we use Extended Euclid 
Algorithm to compute ri , r2 . 

Finally we obtain 



53 = 91 • ^1 + 92 • r2 - (X2 - l)Xf + (5X| - 15^2 + lO)^^ - (8X| - 24^2 + 16)Xi + 4X| - I2X2 + 8 



which vanishes on H and has X'lX2 as its leading term. 

In the same way, we can get gi — {X2 — 1)'^{X2 — 2) 
34 = Xf + 6(X| - 2X2)Xf - 13(X| - 2X2)Xl + 12(X| - 2X2)^1 - 4(X. 



for Xl 



92 



(X2 - l)^Xf for Xi^Xl and 
2X2) for xf. In fact we need 
to compute gi, g2, 93 and 54 in turn according to the lexicographic order because we need reduce 32 by 
gi, reduce 53 by 32 and gi, and reduce 54 by gi, 32 and 33. 

The reduced polynomial set can be proved in Section 6 to be the reduced Grobner basis of the inter- 
section of two ideals which is exactly the vanishing ideal over H, and D is the quotient basis. 
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3. Notions 



First, we define the following mappings. 



proj : Ng — 


k 




_ (rfi,... 


,dn) 


> dn- 


proj : Nq — 






(rfi,... 


, d„) — 


(di, . . . , d„_i). 


embedc ■ Np 






(di,... 


,C?n-l) 


— > (di,. ■ ■ ,(i„_j_^) 



Let D C iVj^*, and naturally we define proj{D) — {proj{d)\d e and embedc{D ) = {embedc{d)\d £ 
D } where D C iVg'"^. In fact we can apply these mappings to any set O C fc" or any matrix of n columns, 
because there is no danger of confusion. For example, let M be a matrix of n columns, and proj{M) is a 
matrix of n — 1 columns with the first n — 1 columns of M reserved and the last one eliminated. 

The embedc mapping embeds an n — 1 dimensional lower set into the n dimensional space. When the 
embedc operation parameter c is zero, we can get an n dimensional lower set by mapping each element 
d = (di, . . . , dn-i) to d = (c?i, . . . , 0) as showed below. 



embedo 




Fig. 6: Embed the lower set in 2-D space into 3-D space with parameter c — 0. 



Blank circles represent the elements of the limiting sets. Note that after the embedc mapping, there is 
one more blank circle. In this case, the limiting set is always increased by one element (0, . . . , 0, 1). 

In the case the embedc operation parameter c is not zero, it is obvious that what we got is not a lower 
set any more. But there is another intuitive fact we should realize. 

Theorem 1: DQ,Di,...,Dk are n — 1 dimensional lower sets, and Dq 3 -Di 3 ... 3 D^. Let 
Di = embedi{Di), i = 0, . . . , fc. Then D — UiLo dimensional lower set, and E{D) C C where 

C = \jtQembed,iE{Di))[j{{0,...,0,k+l)}. ^ ^ ^ ^ 

Proof: First to prove I? is a lower set. Vd G D, let i = proj{d), then d E Di i.e. proj{d) S proj{Di) — 
Di. Because Di is a lower set, hence for j — 1, . . . , n — 1, if dj 7^ 0, then proj{d) — proj(ej) £ Di where 
ej = (0, . . . , 0, 1, 0, . . . , 0) with the 1 situated at the j-th position. So d — Cj G Di C D. For j = n, if 
1 = 0, then we are finished. Else there must be d — e„ e Di-i C D. Because if d — e„ ^ Di~i, we have 
proj{d) ^ Di^i. Since we already have proj{d) G Di, this is contradictory to Di C Z?i_i. 

Second, Vd G E{D), proj{d) ^ Di, i = 0, . . . , k. If proj{d) is a zero tuple, then d„ must be fc + 1, 



that is d G C. Else we know d„ < fc + 1. If d^ 7^ 0, j = 1, 



1 , then d 



G embeddJDaJ. 



Then proj{d) ~ proj{ej) G £'d„, that is proj(d) G E[Dd^). Finally with the embedd^ operation we have 
d G embedd„{E{Dd,J) where d„ < fc + 1. So d G C. 



4. Addition of lower sets 

In this section, we define the addition of lower sets which is the same with that in [2] , the following 
paragraph and Fig. 7 are basically excerpted from that paper with a little modification of expression. 

To get a visual impression of what the addition of lower sets dose, look at the example in Fig.7. What 
is depicted there can generalizes to arbitrary lower sets Di and D2 in arbitrary dimension n and can be 
described as follows. Draw a coordinate system of Ng and insert Di. Place a translate of D2 somewhere on 
the X2-axis. The translate has to be sufficiently far out, so that Di and the translate D2 do not intersect. 
Then take the elements of the translate of D2 and drop them down along the X2-axis until they lie on top 
of the elements of Di. The resulting lower set is denoted by Di + D2. 
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Fig. 7: Addition of Di and L'2- 



Intuitively, we define algorithm AOL to realize the addition of lower sets. 

Algorithm AOL: Given two n dimensional lower sets Di,D2, determine another lower set as the 
addition of Di,D2, denoted by := Z3i + £>2- 
[step 1]: D ■= Di; 

[step 2]: If = return D. Else pick a e 1)2,^^2 := £'2 \ {a}- 

[step 2.1]: If {{D[j{a})=^D, add the last coordinate of a with 1. Go to [step 2.1]. Else D := Dy]{a}, 
go to [step 2]. 

Given n dimensional lower sets Di, D2, D^^ the addition we defined satisfies: 

(a) Di+D2 = D2 + D1, 

(b) {Di + D2) +D3^Di + {D2 + Ds), 

(c) Di + 1)2 is a lower set, 

(d) + D2) ^ ^Di + ^D2. 

These are all the same with that in [2]. And the proof can be referred to it. 

As implied in the example of Section 2, when we want to get a polynomial with leading term ^3 
showed in the right part of Fig. 8, we need two polynomials with the leading terms di,d2 which are not 
the elements of the lower sets and have the same degrees of Xi as ^3 and the minimal degrees of X2 as 
showed in the left part of Fig. 8. In other words, di ^ Di, ^2 ^ D2, proj{di) = proj(d2) = projlds), 
proj{di) +proj(d2) = proj^d^). It's easy to understand that these equations hold for the addition of three 
or even more lower sets. 
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Fig.8: proj[di) = proj{d2) = projids), proj{di) +proj{d2) = projids). 



We use algorithm GLT to get the leading terms di and d2 from ^3 respectively. 

Algorithm GLT: Given a G Ng , and an n dimensional lower set D satisfying a ^ D. Determine another 
r = (ri, . . . , r„) £ Wq which satisfies that r ^ D, proj{r) — proj{a) and (ri, . . . , r„_i, r„ — 1) G D, denoted 
by r := GLT{a,D). 

[step 1]: Initialize r such as proj{r) = proj{a) and proj{r) = 0. 

[step 2]: if r ^ D, return r, else r„ := r„ + 1, go to [step 2]. 

Then di = GLT{d3,Di), da = GLT{dj„D2). 

Definition 3: For any / G k[X], view it as an element in k{Xn)[Xi, . . . , Ar„_i] and define LCn{f) to 
be the leading coefficient of / which is an univariate polynomial of Xn. 

Algorithm GLP: Z? is an n dimensional lower set, a G Nq and a ^ D, G :— {f G k[X]; 3 ed G E{D), s.t. 
the leading term of / is ed }, algorithm GLP returns a polynomial p in the ideal (G) whose leading term 
is GLT{a, D). Denoted by p := GLP{a, D, G). 
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[step 1:] c := GLT{a,D). 

[step 2:] Select c e E{D),s.t. c is a factor of c. d := 4. 

[step 3:] p := f^> ■ d where f^> is an element of G whose leading term is c . 

Remark 1: LCn{f^') = LCn{p) in [step 3]. Since c has the minimal degree of X„ according to 
algorithm GLT, there exists no element c G E{D) which is a factor of c satisfying proj{c ) < proj{c). 
Hence monomial d in the algorithm does not conclude the variable X„. 

5. Associate a lower set D{H) to a set of points H with multiplicity structures 

For any given set of n dimensional points H with multiplicity structures, we can construct an n 
dimensional lower set D{H) by induction. 

Univariate case: H = {{pi, Di), . . . , (pt, Dt)}, then the lower set is D{H) = {0, 1, . . . , tt-^i}- 
To pass from n — 1 to n (n > 2), we first solve a Special case. 

Special case: H = {(pi, Di) , .... (pj, -Dj)} is a sot of n dimensional points with multiplicity structures 
where all the points share the same X„ coordinates. Write H in matrix form as {V ,7)) and all the 
entries in tlic last column of matrix V have the same values. Classify the row vectors of {V^T)) to get 
{(Po:2?o), ■ • • , {'Pw,'Dw)} according to the values of the entries in the last column of matrix V and we 
guarantee the corresponding relationship between the row vectors of matrix V and matrix 2? holds in 
l(Pi,T>i) (0 < i < w). All the entries in the last column of are the same i and the entries of the last 
column of Vi stay the same too. Then eliminate the last columns of Vi and Vi to get {proj{Vi),proj{'Di)) 
which represents a set of n — 1 dimensional points with multiplicity structures, by induction we get a lower 
set Di in n — 1 dimensional space. Then we set 

w 

D{H) = IJ embedi{Di). 

i=0 

Next we deal with the General case. 

General Ccise: H = {{pi, Di) , . . . , {pt, Dt)} is a set of n dimensional points with multiplicity struc- 
tures. Split the set of points: H = Hi[J H2 IJ . . . IJ -ff^- The points of Hi are in the same tt- fibre, i.e. they 
have the same X„ coordinates Cj, i = 1, . . . , s,and 7^ Cj,\li,j = 1,. . . ,s,i ^ j. According to the Special 
case, for each i = 1, . . . , s, we can get a lower set D{Hi), then we set 

s 

D{H) = Y,D{Hi). 

We now proof D{H) is a lower set although it is easy to understand as long as the geometric interpre- 
tation involves. Since it is obviously true for Univariate case, induction over dimension would be helpful 
for the proof. 

Proof: Assume D[H) is a lower set for the n—1 dimensional situation and now we prove the conclusion 
for n dimensional situation (n > 2). 

First to prove D{H) of the Special case is a lower set. 

We claim that {proj{Pi),proj{T>i)) represents an n — 1 dimensional set of points with multiplicity 
structures {i = 0, For any D c Ng, define Fa{D) = {d G D\ proj{d) = a}. Let U = {u\u e 

{1, . . .,t},Fi{Du) 7^ 0}. So {proj{T'i),proj{Vi)) can be written in the form of {(proj(p«),PTOj(-P»(-Du)))|w e 
[/}. Apparently proj{Fi{Du)) is an n — 1 dimensional lower set and can be viewed as the multiplicity struc- 
ture of the point proj(jpu)- Hence {proj{Vi),proj(T>i)) is an n—1 dimensional set of points with multiplicity 
structures. 

What's else, we assert proj{Vj) is a sub-matrix oiprojiVi), and proj{'Dj) is a sub-matrix oiproj{'Di), < 
i < j < w. Because of the corresponding relationship between the row vectors in V and V, we need only to 
prove projiT>j) is a sub-matrix o{projiT>i). If it is not true, there exists a row vector g of proj{'Dj) which 
is not a row vector of proj{'Di). That is, there exists h [1 <h < t) such that emhedj{g) is an clement of 
the lower set and embedi{g) is not included in any lower set Da (1 < a < t). However since i < j and 
embedj{g) G Db, embedi{g) must be included in Di,. Hence our assertion is true. 

Since proj{Pj) is a sub-matrix of proj{Vi), and projiVj) is a sub-matrix of proj{'Di),0 < i < j < w. 
According to the assumption of induction and the way we construct D{H), we have Dj D Dj,0 < i < j < w, 
where Di,Dj are both lower sets. Based on the Theorem 1 in Section 3, D{H) — [Jf^Qembedi{Di) is a 
lower set, and E{D{H)) C (J^^ embedi{E{Di)) U{(0, . . . , 0, w + 1)}. 
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Then to prove D{H) of General case is a lower set. Since D{Hi),i = 1, . . . , s are lower sets, and the 
addition of lower sets is also a lower set according to Section 4, D{H) is obviously a lower set. The proof 
is finished. 

6. Associate a set of polynomials poly{H) to D{H) 

For every lower set constructed during the induction procedure showed in the last section, we associate 
a set of polynomials to it. 

We begin with the univariate case as we did in the last section. 

P-univariate case: H = {(pi, Di), . . . , (pt, A)}, and D{H) = {0, 1, . . . , X;*=i ttA}- The set of 
polynomials associated to D{H) is poly{H) = {11^=1 (-'^i "Pi)"^'}- 

Apparently, poly{H) of P-univariate case satisfies the following Assumption. 

Assumption: For any given n — 1 (n > 1) dimensional set of points H with multiplicity structures, 
there are the following conclusions. For any A G E{D{H)), there exists a polynomial fx G A;[X] where 
X = {Xi, . . . , Xn-i) such that 

• The leading term of fx under lexicographic ordering is X^. 

• The exponents of all lower terms of fx lies in D{H). 

• fx vanishes on H. 

• polyiH) ^ {fx\X e E{D{H))}. 

When we construct the set of polynomials poly{H), we should make sure the assumption always holds. 
Now let us consider the n (n > 1) dimensional situation and still begin with the special case. 

P-Special case: Given a set of points with multiplicity structures H = {{pi, Di), . . . , {pt, Dt)} or in 
matrix form {P = {Pij)mxniT^ = {dij)mxn)- AH the given points have the same X^ coordinates, i.e. the 
entries in the last column of V are the same. We compute poly{H) following the steps below. 

[step 1]: c := w = max{din. i = 1, . . . , m}. 

[step 2]: Vi = 0, . . . , w, define SVi as a sub-matrix of V containing all the row vectors whose last coor- 
dinates equal to i. Extract the corresponding row vectors of V to form matrix SVi, and the corresponding 
relationship between the row vectors in V and T) holds for SV i and SVi. 

[step 3]: Vi = 0, . . . , w;, eliminate the last columns of SVi and SVi to get {SVi, SVi) which represents a 
set of points in n— 1 dimensional space with multiplicity structures. According to the induction assumption, 
we have the polynomial set Gi = poly {{SVi, SVi)) associated to the lower set Di = D{{SVi,SVi)). 

[step 4]: D := [j^^Qembedi{Di). Multiply every element of Gj with (X„ — c)* to get Gj. G := 

ur=oG.u{(^n-cr+^}- 

[step 5]: Eliminate the polynomials in G whose leading term is not included in E(D) to get poly{H). 
Theorem 2: The poly{H) got in P-Special case satisfies the Assumption. 

Proof: According to the Section 5, {SVi,SVi) represents an n — 1 dimensional set of points with 

multiplicity structures for i = 0, . . . ,w. And Dj D Di,0 < j < i < w. D is a lower set and E{D) C 
Ur=o embed,{E{bi)) U{(0, ...,Q,w + 1)}. 

For A = (0, . . . , 0, w + 1) e E{D), we have fx = [Xn — c)"'+^. It is easy to check that it satisfies the 
first three terms of the Assumption. 

For any other element ed of E{D), 3k s.t. ed G embedkE{Dk). So let ed be the element in E{Di.) such 
that ed = embedk{ed). We have f~^ vanishes on {SVk,SVk) whose leading term is ed G E{Dk) and the 
lower terms belong to Dk- According to the algorithm fed = (X„ — c)'^ • /g-^ G poly{H) . 

First it is easy to check that the leading term of fed is ed since ed = embedk{ed). 

Second, the lower terms of fed are all in the set S — Ui=o embed j{Dk) because all the lower terms of 
f-^ are in the set Dk- Dq D Di D . . . Dk, so embed j{Dk) C embed j{Dj) {0 < j < k), hence S C D = 
Uj=o embed j{Dj) and the second term of the Assumption is satisfied. 

Third, we are going to prove that fed vanishes on all the functionals defined by {V,V), i.e. all the 
functionals defined by {SVi, SVi) {i — 0, . . . ,w). 

When i ^ k, we write all the functionals defined by {SVi, SVi) in this form: L' ■ g^|x„=c where L' 
is an n — 1 variable functional. Since fed — {Xn — c)^ ■ fed: apparently fed vanishes on these functionals. 

For i = k, denote by L the functionals defined by {SVi, SVi), and /^^ vanishes on L. All the functionals 
de&ned hy {SVk,SVk) can be written in this form: L "• |x„=c where L" G L. Since /erf = (X„ — c)'^-/g-^, 
apparently fed vanishes on these functionals. 

So fed vanish on H, and fed satisfies the first three terms of the Assumption. 

In summary poly{H) satisfies the Assumption, and we finish the proof. 



8 



Remark 2: For fx G poly{H).X e E{D) where poly{H) is the result got in the algorithm above, we 
have the conclusion that LC„(/a) (X„ - c)p™^(^\ 

P-General case: Given a set of points with multiplicity structures H or in matrix form {V = 
{Pij)mxn,'D = {dij)mxn), we are going to get poly(H). 

[step 1]: Write H as H = Hi[j H2[j ■ ■ - {J Hg where -ffj (1 < i < s) is a 7r-fibre {-k : H ^ k such that 
{p = {pi, . . . ,Pn), D) e iJ is mapped to Pn € k) i.e. the points of Hi have the same X„ coordinates Cj, 
i = 1, . . . , s,and Ci Cj,\/i,j = 1, . . . , .s, i ^ j. 

[step 2]: According to the P-Special case, we have = D{Hi),Gi = poly{Hi). Write Hi as {Vi,T>i), 
and define Wi as the maximum value of the elements in the last column of P,. 

[step 3]: D := D[,G := Gi,i := 2. 

[step 4]: If i > s, go to [step 5]. Else 

[step 4.1]: D := D + Df, G := 0. View E{D) as a monomial set MS := E{D). 

[step 4.2]: If ^MS = 0, go to [step 4.7], else select the minimal element of MS under lexicographic 
ordering, denoted by LT. MS := MS \ {LT}. 
[step 4.3]: 

/i ■.^GLP{LT,D,G),h ■.= GLP{LT,D\,Gi). 
Vk ■■=proj{gk), where gk := GLT{LT, D',^), k = l,...,i. 

[step 4.4]: 

i-l 

fe=i 

i-l 

ppl := (X„ - Cir^+^--- PP2 := n(^n - c,)-'=+l-''^ 

fe=i 

[step 4.5]: Use Extended Euclidean Algorithm to compute ri and r2 s.t. ri ■ ppi + r2 ■ pp2 = 1- 
[step 4.6]: f := ri ■ qi + r2 ■ q2- Reduce / with the elements in G to get f ; G := G\J{f }. Go to [step 
4.2]. 

[step 4.7]: G := G. i ■.= i + l. Go to [step 4]. 
[step 5]: poly{H) := G. 

Theorem 3: The poly{H) got in p-General case satisfies the Assumption, 
proof: Wc need only to prove the situation that s > 2 in [step 1]. 

For i = 2, D = Di + ^Dj- e E{D), v := proj{ed) and Xq According to Section 4, we 

have V = vi + V2- Based on the Remark 1 and Remark 2, /i and /2 can be written as polynomials of 
k{Xn)[Xi, . . . , Xn-i] : fi = Xo ■ {Xn — Ci)"! + the rest and f2 = Xq- (X„ — 02)"^ + the rest and none of 
the monomials in the rest is greater than or equal to Xq. Because /i and (X„ — ci)"'^^^ vanish on Hi, f2 
and {Xn - C2)"'^+^ vanish on H2, we know that qi = fi ■ (X„ - 02)""^+^ and q2 = /2 • (X„ - ci)"'i+i both 
vanish on Hi\jH2. Then / vanishes on Hi [J H2 where f = ri ■ qi + r2 ■ q2- 

/ = Xo • (X„ - ciY^ ■ {Xn - C2r-{ri ■ {Xn " 02^^+'""=^ + ^2 ' - Ci)-i+i-''0 + t/^G rest 
= Xo ■ {Xn - CiY' ■ {Xn - 02)"= (ri • ppl + r2 ■ pp2) + the rest 
= Xo ■ {Xn - cifi • {Xn - C2y^ + the rest 

None monomial in the rest is greater than or equal to Xq , so the leading term of / is apparently 
Xq ■ Xn which is equal to ed. Moreover we naturally have the following Proposition 1 for i = 2. 

Proposition 1: For every polynomial / we get in the algorithm, LC„(/) = 11^=1 (^n ~ ^jT^ ■ 

When i > 2, assume the Proposition 1 holds for i — 1. V ed e E{D), v := proj{ed) and Xq := 
According to Section 4, we have v = vi + . . . + Vi. Based on the Proposition 1, Remark 1 and Remark 
2, fi and /2 can be written as polynomials of A;(X„)[Xi, : /i = Xq ■ '[f'~}^{Xn — cj)'"^ + the rest 
and f2 = Xo ■ {Xn — CiY* + the rest and none of the monomials in the rest is greater than or equal to 
Xq. Because /i and ]~[j=i(^n — Cj)'^''^^ vanish on Uj=i (X^ — Ci)™'+^ vanish on we know 

that qi = fi ■ {Xn — Ci)"'*+^ and ^2 = /2 • 11}=! ~ both vanish on U}=i ^j- Then / vanishes 

oil Uj=i where f = ri ■ qi + r2 ■ q2- 

i i — 1 

f = Xo- l[{Xn - cjP{ri ■ {Xn - + r2 • l[{Xn - cjr^+^-^^) + the rest 

= Xo ■ ]Tj=i{Xn - c^Y' (r-1 • ppl + r2 • pp2) + the rest 
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= ^0 • IVj^ii^n - Cjp + the rest 

None monomial in the rest is greater than or equal to Xq and the leading term of / is apparently 
Xq ■ X^ which is equal to ed. Hence the Proposition 1 holds for arbitrary i. 

Therefore we have proved that for any element ed G E{D), fed := f vanishes on H and the leading term 
is ed. In the algorithm, we compute fed in turn according to the lexicographic ordering of the elements of 
E{D). Once we get a polynomial, we use the polynomials we got previously to reduce it ([step 4.6]). Now 
to prove the lower terms of / are all in D after such a reduction operation. 

Let D be a lower set, a be a monomial, define L{a,D) = {b £ No;6 ^ a,b £ D}. Given any d ^ D, 
there exist only two situations: d G E{D) or ^ E{D) but 3d G E{D), s.t. is a factor of d. Of course 
d' -< d. 

The very first vanishing polynomial we got in the algorithm is an univariate polynomial of X„ with 
leading term being T. It is easy to check it's lower terms are in D. Sincx; the polynomial is a vanishing 
polynomial, we can say that T can be represented as the linear combination of the elements of L(T, D). 

Since T is the first element which is not in D under lexicographic ordering. We assume that there 
exists such a monomial M ^ D{M >- T) that Vm -< M{m ^ D), m can be represented as the linear 
combination of the elements of L(m, D). Now to prove M could be represented as the linear combination 
of the elements of L{M, D). 

If M € E{D), then the algorithm provides us a vanishing polynomial whose leading term is M i.e. 
that M can be represented as the combination of the terms which are all smaller than M. According to 
the assumption, for any lower term m {m ^ D) of the polynomial, m can be represented as the linear 
combination of the elements of L{m,D), then M could be represented as the linear combination of the 
elements of L{M,D). 

li M ^ E{D), there exists rf' e E{D) s.t. M = M' ■ d! . Since d! -< M, according to the assumption, 
we can substitute d with the linear combination of the elements of L{d ,D). Since all the elements in 
L{d , D) are smaller than d , then M could be represented as the combination of elements which are all 
smaller than M. Then for the same reason described in the last paragraph, M could be represented as the 
linear combination of the elements of L{M, D). 

Therefore specially for any ed e E{D), all the lower terms of the polynomial fed we got in the algorithm 
after the reduce operation are in D, and the proof is done. 

Theorem 4: Given a set of points H with multiplicity structures, poly{H) is the reduced Grobner 
basis of the vanishing ideal I{H) and D{H) is the quotient basis under lexicographic ordering. 

Proof: Let m be the number of functionals defined by H and then m = dim{k[X]/ I{H)). Denote by J 
the ideal generated by poly{H). According to the Assumption, poly{H) C I{H). So dim{k[X]/I{H)) < 
dim{k[X]/J). Let C be the leading terms of polynomials in J under lexicographic ordering, then C 3 
U/3eis(D(i/))(/3 + No) where the latter union is equal to Ng-D(ff). Then we can get C' = N^-C C D{H). 
Because k[X]/J is isomorphic as a A;- vector space to the fc-span of C , here C is viewed as a monomial 
set. So dim{k[X]/J) < ^D{H) = m. Hence we have 

m = dim{k[X]/I{H)) < dim{k[X]/J) < m. 

Therefore J = I{H), where J = {poly{H)). Hence apparently poly{H) is exactly the reduced Grobner 
basis of the vanishing ideal under lexicographic ordering, and D{H) is the quotient basis. 

7. Intersection of ideals and some applications 

Some steeps of our algorithm actually do the work of computing the intersection of two ideals, but 
we note that the information of the zeros of the ideals is necessary there (see [step 4.1] - [step 4.7] of 
p-General case in Section 6). We now bring up a new algorithm to compute the intersection of two 
ideals which does not require the information of the zeros of the ideals. 

Lemma 1: G is the reduced Grobner basis of some n- variable polynomial ideal under lexicographic 
ordering with Xi >- X2 y . . . >- X„. Define po(G) as the univariate; polynomial of X„ in G. View g G G 
as polynomial of K{Xn)[Xi, . . . ,X„_i] and define LGn{g) to be the leading coefficient of g which is an 
univariate polynomial of X„ and we have the conclusion that LGn{g) is always a factor of po{G). 

Proof: In fact Proposition 1 in Section 6 holds for any given reduced Grobner basis under lexico- 
graphic ordering since it is unique and can be constructed in the way our algorithm offers. According to 
the proposition, V/ € G, LCn{g) = Uj^A^n - Cjf^ and < w, + 1. po(G) = nj=i(^« - Cj)'"^+^ Hence 
the proof is done. 

Based on Proposition 1 and Lemma 1, we give the algorithm Intersection to compute the in- 
tersection of two ideals /i and I2 which are represented by the lexicographic ordering reduced Grobner 
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bases Gi and G2 and the greatest common divisor of po(Gi) and ^0(^2) equals to 1. Denote by Q{G) the 
quotient basis where G is the reduced Grobner basis. Algorithm GP is a sub-algorithm called in algorithm 
Intersection. 

Algorithm GP: G is a reduced Grobner basis, for any given monomial LT which is not in Q{G), we 
get a polynomial p in (G) whose leading term is a factor of LT: the Xi,. . . ,Xn-i components of the 
leading term arc the same with that of LT and the X„ component has the lowest degree. Denoted by 

p:= GP{LT, G). 

[step 1:] G := {5 G G| the leading monomial of 5 is a factor of LT }. 

[step 2:] G := {g & G \'^g e G , s.t. the degree of X„ of the leading monomial of g is lower than 
that of g }. 

[step 3:] Select one element of G and multiply it by a monomial of Xi, . . . to get p whose 

leading monomial is LT. 

Algorithm Intersection: Gi and G2 are the reduced Grobner bases of two different ideals satisfying 
that GCD{po{Gi),po{G2)) = 1. Return the reduced Grobner basis of the intersection of these two ideals, 
denoted by G := Intersection{Gi,G2)- 

[step 1:] D := Q(Gi) + Q{G2). View E{D) as a monomial set. G := 0. 

[step 2:] If E{D) = 0, the algorithm is done. Else select the minimal element of E{D), denoted by T. 
E{D) := EiD)/{T}. 
[step 3:] 



[step 4:] 



= GP(r,Gi), /2 :=GP(r,G2). 
qi ■■= fi ■P2, 92 := /2 - Pi- 



. Po(G2) ^ PoiGi 



LCM2)' LCnih)' 

[step 5:] Use Extended Euclidean Algorithm to find ri,r2 s.t. 

ri-ti+r2-t2 = 1. 

[step 6:] / (ji • ri + (j2 • ?'2- Reduce / with G to get / , and G := G[J{f }. Go to [Step 2]. 

Because the algorithm is essentially the same with [step 4.1] - [step 4.7] of p-General case in Section 
6, here we don't give the proof. 

The Proposition 1 and Lemma 1 reveal important property of the reduced Grobner basis under 
lexicographic ordering. If a set of polynomials does not have this property, it is surely not a reduced 
Grobner basis. 

It is well-known that the Grobner basis of an ideal under lexicographic ordering holds good algebraic 
structures and hence is convenient to use for polynomial system solving. To compute the zeros of an zero 
dimensional ideal with the reduced Grobner basis G, we need first compute the roots of po{G). Since 
LCn{g) {g 7^ Po{G), g G G) is a factor of po(G), compute the roots of LCn{g) which has a smaller degree 
would be helpful for saving the computation cost. 



8. Conclusion 

Based on the algorithm Intersection in Section 7, the algorithm of p-General case in Section 6 can 
be simplified. The last sentence in [step 2] can be deleted and we can replace [step 4.3] and [step 4.4] by: 
[step 4.3]: 

/i := GLP{LT,D,G),f2 := GLP{LT, D^Gi). 

[step 4.4]: 

qi ■■= fi ■ PoiGi); 92 := /2 • Po{G). 



~ LCnihY LCMiY 

During the induction of the algorithm in Section 6, we can record the leading coefficients for later 
use to save the computation cost and the computation cost is mainly on the Extended Euclid Algorithm. 
However the advantage of our algorithm is not fast computation, after all it depends on how many times 
we need to use the Extended Euclid Algorithm. 

Our algorithm has an explicit geometric interpretation which reveals the essential connection between 
the relative position of the points with multiplicity structures and the quotient basis of the vanishing ideal. 
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The algorithm offers us a new perspective of view to look into the reduced Grobner basis which can help 
us understand the problem better. Lemma 1 and the algorithm to compute the intersection of two ideals 
are the direct byproducts of our algorithm. 

Since we finished the paper [1] previously which gives an algorithm to get the minimal monomial basis 
of Birkhoff interpolation problem with little computation cost, we have always believed that the algorithm 
could be interpreted in a more geometric way and the proof should be more beautiful and much easier to 
understand. The proof in [1] is so complicated that we ourselves don't like it. And it would be great if 
we can get the interpolation polynomial with little computation cost instead solving the linear equations 
since the minimal monomial basis can already be got in a simple way. That's why we began to study 
the vanishing ideal of the set of points with multiplicity structures which is essentially a special case of 
Birkhoff interpolation problem. 

I still remember the moment when I first read the paper [2] written by Mathias Lederer in which the 
quotient basis and Grobner basis can be got in a geometric way. I told myself that this was just what we 
wanted. Paper [2] concentrates on the the vanishing ideal of the set of points with no multiplicity structures 
in afBne space. Although whether or not the points are with multiplicity structures matters much, the 
paper really inspired us a lot. Our algorithm also uses induction over variables and the definition of addition 
of lower sets is essentially the same with that in paper [2] . However during the induction procedure, we 
have to consider p-Special case and p-General case. This consideration, on one hand, clearly indicates 
the geometric meaning of the multiplicity structures of points, on the other hand, means a lot for our 
capacity of applying the algorithm of intersection of two ideals. In paper [2], the author uses Lagrange 
interpolation method to get the vanishing polynomial of all points from the polynomials vanishing on 
subsets of the points. However the Lagrange interpolation method could just not work for our problem 
because the points are with multiplicity structures. In this paper, we creatively use the Extended Euclidean 
Algorithm. Thanks goes to paper [2] and the author, after we solved the problem of the vanishing ideal of 
the set of points with multiplicity structures, we will move on to the Birkhoff problem. 
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